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Abstract. In this paper we show that the h-p spectral element method developed in 
1 3 8 9 1 applies to elliptic problems in curvilinear polygons with mixed Neumann and 
Dirichlet boundary conditions provided that the Babuska-Brezzi inf-sup conditions are 
satisfied. We establish basic stability estimates for a non-conforming h-p spectral ele- 
ment method which allows for simultaneous mesh refinement and variable polynomial 
degree. The spectral element functions are non-conforming if the boundary conditions 
are Dirichlet. For problems with mixed boundary conditions they are continuous only 
at the vertices of the elements. We obtain a stability estimate when the spectral element 
functions vanish at the vertices of the elements, which is needed for parallelizing the 
numerical scheme. Finally, we indicate how the mesh refinement strategy and choice of 
polynomial degree depends on the regularity of the coefficients of the differential oper- 
ator, smoothness of the sides of the polygon and the regularity of the data to obtain the 
maximum accuracy achievable. 
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1. Introduction 

In this paper we generalize all the results we have obtained in 1 3 1 and seek a numerical 
solution to an elliptic boundary value problem where the differential operator satisfies the 
Babuska-Brezzi inf-sup conditions. We solve the boundary value problem on a curvilin- 
ear polygon whose sides are piecewise analytic (smooth) and we assume the boundary 
conditions are of mixed Neumann and Dirichlet type as in I1I2I5I . 

We now briefly describe the contents of this paper. In §2 we discuss function spaces 
and obtain differentiability estimates for the solution in modified polar coordinates in a 
sectoral neighbourhood of the vertices. Here we examine two cases viz. when the coeffi- 
cients of the differential operator, sides of the polygon and the data are analytic and when 
they have finite regularity. 

In §3 we obtain a stability theorem for a non-conforming spectral element represen- 
tation of the solution for problems with mixed boundary conditions. We let the spectral 
element functions to be polynomials of variable degree, where the degree of all these poly- 
nomials is bounded by W, and let M denote the number of elements or layers in a sectoral 
neighbourhood of each of the vertices in the radial direction as shown in figure^ We then 
define a quadratic form y M -W which measures the sum of squares of a weighted squared 
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norm of the partial differential equation and fractional Sobolev norms of the boundary 
conditions and a term which measures the jumps in the function and its derivatives at 
inter-element boundaries in appropriate Sobolev norms. In each of the sectoral neighbour- 
hoods of the corners we use modified polar coordinates and a global coordinate system in 
the remaining part of the domain. We prove that the sum of the squares of the H 2 norms 
of the spectral element functions is bounded by the quadratic form y M < w multiplied by 
a factor which grows logarithmically in W for problems with Dirichlet boundary condi- 
tions. For problems with mixed boundary conditions this factor can grow as M 4 , provided 
W is not too large, and thus the method displays algebraic instability. 

We choose as our approximate solution the unique spectral element function which 
minimizes a functional v 3 ^^ closely related to the quadratic form y M ' W as defined in 
I3I8I9I . In case the solution is analytic, we choose M proportional to W, and show that 
r 14 ^ decays exponentially in M. Now the error is bounded by r M W multiplied by a factor 
which grows at most algebraically in M. Hence the order of convergence remains expo- 
nential. If the solution has finite regularity then we choose M proportional to InW and 
show that r 14 ^ decays algebraically in W. Now the error is bounded by r™' W multiplied 
by a factor which grows polylogarithmically in W and hence the error decays algebraically 
in W. 

We now come to the aspect of parallelization of the numerical scheme. For problems 
with Dirichlet boundary conditions the spectral element functions are non-conforming 
and we can use the stability theorem to parallelize the scheme in an optimal manner. 
It should be noted that the method is assymptotically faster then the h-p finite element 
method. For problems with mixed boundary conditions we cannot use this stability the- 
orem to parallelize our method since the factor in the stability estimate can grow as M 4 . 
To get around this problem we make the spectral element functions continuous at the ver- 
tices of the elements only. We then prove a stability theorem for mixed problems when 
the spectral element functions vanish at the vertices of their elements. The values of the 
spectral element functions at the vertices of their elements constitute the set of common 
boundary values we have to solve for. It should be noted that the cardinality of the set 
of common boundary values is much smaller than for finite element methods where the 
functions have to be continuous along the edges of the elements. Since the cardinality of 
the set of common boundary values is small we can construct an accurate approximation 
to the Schur complement matrix from its definition. As a result the method is faster than 
the standard h-p finite element method JSJ . 

2. Function spaces and differentiability estimates 

Let f2 be a curvilinear polygon with vertices Aj,A2, . . . ,A p and corresponding sides 
r^T^, . . . ,T p where T,- joins the points A,_i and A, . We shall assume that the sides T, 
are analytic (smooth) arcs, i.e. 

r ; = {(<?,-(§), y/^))|£e7=[-U].} 

with <p,-(|) and V / i(^) being analytic (smooth) functions on 7 and |<p/(<i;)| 2 + | V / /(^)| 2 > 
a > 0. By r, we mean the open arc, i.e. the image of / = (—1,1). Let the angle subtended 
at Aj be 03j. We shall denote the boundary d£l of D. by T. Further let T = \JT^\ 
r[01 = U, G @-r ; , rW = U; e ,/r T; where S> is a subset of the set {i \ i = 1, . . . ,p} and jY = 
{i | / = 1 , . . . ,p} \ *3. Let x denote the vector x = (x\ ,X2). 
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Let £ be a strongly elliptic operator 

2 2 

£(«) = 52 feW^k + H b r (x)u Xr +c(x)u, (2.1) 

r,s = 1 r= 1 

where fl S) r(x) = a^sM, b r (x), c r (x) are analytic (smooth) functions on Q. and for any 
(|i , I2) S K and any xeH, 



2 

I 

r,.v=l 



E «^ s >Mofe 2 + ^2 2 ) (2-2) 



with /io > 0. Moreover let the bilinear form induced by the operator £ satisfy the inf-sup 
conditions. 

In this paper we shall consider the boundary value problem 

Zu = f on £1, 

u = gW onrM, 

M)r^ onrl '' 

where j denotes the usual conormal derivative which we shall now define. Let A 
denote the 2x2 matrix whose entries are given by 

A rtS (x) = a r . s (x) 

for r,s = 1,2. Let N = (A^A^) denote the outward normal to the curve T, for i e jV , 
Then (§]fj is defined as follows: 

— x = £ A^ — . (2.4) 

We shall assume that the given data / is analytic (smooth) on i2 and is analytic 
(smooth) on every closed arc r, and is continuous on V™ . 

We need to state our regularity estimates in terms of local variables which are defined on 
a geometrical mesh imposed on Q. as in §5 of [2 1. We first divide Q. into subdomains. Thus 
we divide Q. into p subdomains S ,S P , where S' denotes a domain which contains the 
vertex A' and no other, and on each S' we define a geometrical mesh. Let & k = •,_/ = 
l,...,J k ,i= l,...,4j} be a partition of S k and let 6 = U P k=l & k . Here J k =M + 0{\) 
and 4 , / < / for all and 7", where / is a constant. As has been stated earlier M denotes the 
number of elements or layers in a sectoral neighbourhood of each of the vertices in the 
radial direction. 

We now put some restrictions on 6. Let (rj, 9^) denote polar coordinates with center at 
Afc. Let Ik = lnr^.. We choose p so that the curvilinear sector Q* with sides and r^ + i, 
center at A^ and radius p satisfies 

U 
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r, a 




Figure 1. Geometric mesh with M layers in the radial direction in the curvilinear 
domain. 



£2 may be represented as 



£l k = {(x,y)eQ.:0<r k <p}. 

The geometrical mesh we have imposed on £2 is as shown in figure^ 

Let yfj : j, 1 < / < 4 be the side of the quadrilateral Q. k j 6 &. Then we assume that 



(2.5) 



i ■ 



0<% < 1,/ = 1,3 



0< tj < 1,1 =2,4 



and that for some C > 1 and L > 1 independent of i,j,k and / 



d' 



d' 



< CL't\,t = 1,2,. 



(2.6a) 
(2.6b) 

(2.7) 



We shall also examine the case when they are smooth. Some of the elements may be 
triangles too |7 1. We shall place further restrictions on the geometric mesh we impose on 
D. k later. 




Figure 2. Curvilinear sectors. 



Stability estimates for h-p spectral element methods 



399 



Let (r k , Ok) be polar coordinates with center at A k . Then £t k is the open set bounded 
by the curvilinear arcs T kl T k +\ and a portion of the circle r k = p. We subdivide Q. k into 
curvilinear rectangles by drawing M circular arcs r k = oj' = p/4 J ,j = 2, . . . ,M + 1 , 
where \i k < 1 and 4 — 1 analytic curves . . . , C/, whose exact form we shall prescribe 
in what follows. We define of = 0. Thus 4 ; = 4 for j < M; in fact, we shall let 4,/ = 4 
for /' < M + 1 . Moreover 4 i < I for all A;, j where / is a fixed constant. Let 

r k+J = {(r k ,0 k )\6 k = fj(r k ),0 <r k < p}, 

j = 0, 1 in a neighbourhood of Then the mapping 

r k = Pk, e k = {w ^_ wk) Kfe - Vf)fi(pk) - (<h - v£)/o (P*)L (2.8) 

where ff is analytic in r k for j = 0, 1, maps locally the cone 

{(p*, fe) : < p k < a, y/f < # < y/f} 

onto a set containing Q. as in §3 of 0. The functions /} satisfy f£(0) = y/f, f\ (0) = y/f 

and (/y) (0) = for j = 0,1. It is easy to see that the mapping defined in (12.81 has two 
bounded derivatives in a neighbourhood of the origin which contains the closure of the 
open set 

6 A = {(At, fc) : < p k < p , y/f < t < y/f}. 

We choose the 4-i curves C2, . . . ,C/ t as 
C i :(j> k (r k ,d k ) = \rf 

for / = 2, . . . ,4. Here y/f = y/f < y/* < • • • < y/f +1 = y/*. Let A y/f = y/f +1 - y/f. Then we 
choose {y/f },jt so that 

max(Ay/f) < A(min(Ay/f)) (2.9) 

i,k i,k 

for some constant A. We need another set of local variables (x k , 0^) in a neighbourhood 
of Q. k where T k — \nr k . In addition we need one final set of local variables (y k , <l>k) in the 
cone 

{{pk, h)-0<p k <p,Wi<h< Yu}, 

where V k = \np k . Let 5* = {(r k , k ) : < r k < jj.} n i2. Then the image 5^ in (v*, 0*) 
variables of S 1 ^ is given by 

% = {(V*,^) : — < v, < ln M , y/f < <j> k < y/f}. 

Now the relationship between the variables (z k , 9 k ) and (v^,^) is given by (T k ,9 k ) = 
M k (v k ,ij> k ), viz. 

T* = V fe , 

0A = (^T^ Kfe - Vf)/f (e n ') - (fc - vf)/oV*)l- (2.10) 
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Hence it is easy to see that J M k (Vfc, <pk), the Jacobian of the above transformation, satisfies 
Ci < \J Mk {v k ^ k )\ < C 2 for all {v k ,<j> k ) G S*, for all < /i < p. 

We should mention here that it is not necessary to choose the system of curves we have 
chosen to impose a geometric mesh on S*, . However it is necessary to choose the curve 
r k = p as the boundary of £2* and no other, as will become apparent in what follows. Any 
other additional set of analytic curves which imposes a geometrical mesh on S k ^ would 
do equally well. However the set of curves we have chosen is, in some sense, the most 
natural as the image £2? • of a curvilinear rectangle £2? ■ for j > 2 in (v^, <p k ) variables is 
given by a rectangle with straight lines for sides and for j = 1 is a semi-infinite strip with 
straight lines for sides. 

We now state the differentiability estimates for the solution u of (12. 3> which will be 
needed in this paper. 

PROPOSITION 2.1. 

Consider the case when the coefficients of the differential operator are analytic on £2 and 
the sides of the curvilinear polygon are analytic. Moreover let the geometric mesh satisfy 
&2.7\ . Let the data f be analytic on £2 and let be analytic on every closed arc Ti ,for 
I = 0,1, and let g^ be continuous on r^L Let Uffv^fa) = u(Vk,<Pk)for (V£,<fe) £ ft?/ 
for j < M and a k = u(Ai l ). Now there is an analytic mapping ■ : SarrowQ. k - for j > M 
given by Mf = {X k j (E,,7]),Y k j (E,,r\)). Here S is the unit square. Let Tj) = 

w(Xf ■ , 17 ) , Yf, (§ , TJ )). r/ie« we can i/zow as in ll3l8l fto 

||^(v,,0,)-«,||^ < (CmUrtf-XH"-* 2 ))* (2.11a) 

/or 1 <j<M,k= l,...,p,l < i < Ik and 

\\U^,f 1 )\\ 2 m<s <(Cmld m ) 2 (2.11b) 

for M<j<Jk, l<i< Ikj, 1 < k < /?. i/ere C,t/ anc/ j3,t are constants and < jQ^ < 1 
/or 1 < Jt < p. 

We next consider the case when the data has finite regularity. To state the differentia- 
bility results in this case we shall need to use the space Hp (£2) with k > I defined in 
(T). We now cite Remark 3 after Theorem 2.1 of Q. Let F] eC m+2 (I) for j = l,...,p 

3 3 

and let the coefficients of the differential operator eC'"(£2). Let eH? +2 ' 2 (T^), 
g [i] g ^"+3.5 (p[i]) and f g H%'°(£1). Then there exists a constant tf m such that 



l 



Nl^o, <K m ll/ll <(£J) + £||g^|| 3_,3_ ; ^ r;ix ] • (2.12) 



PROPOSITION 2.2. 

Consider the case when the differential operator and data satisfy the conditions stated 
above. We assume moreover that the curves (j> k j t and y/fy t defined in J2.6a> . (I2.6bi satisfy 

Ui,j,l\\m+2,-,h II Vf^ll*fW ^ E "'+ 2 
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where E m+ 2 is a constant independent of i,j,k and I. Let UfAv^,^) = u(Vk,<j>k) for 
(Vk, §k) £ "f jf or J < M and a/, = M(A,t). Now there is a smooth mapping Mf ■ : SarrowQ^ ■ 
/or j > M given by Mf M,T]) = (X^.(^,Tj),yA(^,Tj)). i/ere 5 M f/ie wnfr square. Let 
JJ*. (| , tj) = 17 (| , TJ ), (§ , TJ )). 77ien itfihg d2~L2l we can s/zow f/iaf 



/or 1 < y < M, = 1 , . . . , p, 1 < i < Ik and 

\\utj(^ri)t + 2,s<K m+2 (2.13b) 
for M < j < J k , 1 < i < 4j, 1 < < p. Here K m+ 2 denotes a constant. 

3. Stability estimates 

3.1 Preliminaries 

Let 

2 2 
£« = - Y i a r,s( X ) U Xs)x r + Y b r (x)u Xr +c(x)u (3.1) 
r,s=l r=\ 

be a strongly elliptic operator which satisfies the inf-sup conditions. Hence there exists a 
positive constant /Xo > such that 

2 

r,s=l 

for all x £ Q.. 

Let # = Hq(Q.) where w e ffj(^) ^ w e H 1 ^) and trace (w)| r|0 ] = 0. Consider the 
bilinear form B(m, v) defined on H x i/ as follows: 

B(u,v)= V a rs (x)u x s v Xr + V b r (x)u Xr v + cuv )dx. (3.2) 

Then B(u,v) is a continuous mapping from // x HarrowK. and there exists a constant Ci 
such that 

|B(«,y)|<C 1 ||«|| // i ( a)lkll/ / i(n) (3-3) 
for all m , v € Hq (£1). Moreover we assume that the inf-sup conditions |7 1 

B(u,v) 

inf sup — >C 2 >0, (3.4a) 



and 



supB(«, v) > for every ^ v e H (3.4b) 

ueH 
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hold. Then for every continuous linear functional F(v) defined on Hq (Q.) there exists 
unique uq e Hq(Q.) such that B(uq,v) = F(y) for all v G Hq(Q). Moreover, the a priori 
estimate 



\ u o\\hUq.) ^ 



1 



L " 2 O^vi 



sup 



l^(v)| 



(3.5) 



holds. 



Now consider the following mixed boundary value problem 
£w = f in £2, 



y Q u = M| r[0 ] 



and 



7i" 



(9m \ 



(3.6a) 
(3.6b) 



(3.6c) 



Here the conormal derivative y x u is defined as follows. Let T,- C rW and let T and AT 
denote the unit tangent vector and unit outward normal at a point P on T, which we 
traverse in the clockwise direction. Let T = (T\, T^f and N = {N\,N2f . Then 



7i"lr, 



du 
~dN 



2 

r, r,s=l 



du , _, 
a r <—— =NAV x u. 
ox s 



(3.7a) 



In the same way we define the cotangential derivative 

du 

' dx s 



du 
df 



= £ T r a rjS p-=T l ,\V rll . 

Ti r,s=l 



(3.7b) 



and the tangential vector 



du 
df 



(3.7c) 



We now consider the spectral elements which are not contained in the sectoral neighbour- 
hoods of the vertices £l k for k = 1 , . . . , p. Now Q, k j C £l k for 1 < i < I k j and 1 < j < M. 
Let 

O p+1 = {Q. k j, 1 < k < p,M < j < J k , l<i< I kj }. 
Once more J k = M + O ( 1 ) . We shall relabel the elements of O p+ 1 and write 



0" +1 = {£lf +1 ,l <1<L}. 



We shall now introduce some notation so that the reader may proceed directly to the 
stability theorem l3~2l and examine the proof later as it is quite involved. 
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Consider the domain Q.f + . Then there is a mapping Mf + from the master square 
S = (0,1) x (0,1) to £lf + . Let jf + (§,T?) denote the Jacobian of the transformation 
Mf +1 . We let 

w w 

«f +1 (l,r/) = L EMV. 

y=o i=o 

We choose the spectral element functions {w* ;{Vki§k)}i,i,k f° r 1 < i < 7t, 1 < J < Af 
and 1 < k < p to be polynomials of the form 

Wj Wj 

for /' 7^ 1. Here 1 < < W. If 7=1 we choose u k n (Vk^k) — 8k where is a constant for 
1 < i < h- Let n M ' W denote the space of polynomials {{wf\(Vfc, <\>k)}i,j,ki { u f +l (4 , ? ?)}z}- 

Remark 1. We shall always choose M = 0(W). In case the conditions of Proposition l2.ll 
are satisfied so that u is analytic we choose W = M. Once we have obtained the numerical 
solution we can define a correction to it so that the corrected solution is conforming and 
converges to the actual solution exponentially in M in the H 1 (Q.) norm [8 9|. Thus the 
error in the H l (Q.) norm is bounded by Cer bM where C and b are constants. In case 
u e H 1 ? ' (£2) we would choose M proportional to mlnW. Once more we can define 
a corrected version of the solution so that it is conforming and converges to the actual 
solution in the H l (£2) norm and the error is bounded by C(lnW) 3 W~ m+1 . Hence for the 
method to converge we must have m>2. 

The stability theorem 13.21 holds provided the coefficients of the differential operator 
G C 3 (£2) and the curves 0* ■ z , i//f • , defined by ( I2.6ab . J2.6bl > satisfy 

llC/lk^l!^;,/H3W<^3, 

where K3 is a constant independent of i,j,k and I. In this paper however we prove The- 
orem !3 .21 assuming that the coefficients of the differential operator are analytic on £2 and 
the curves fof . v/f . , defined in ( I2.6ab . J2.6bi are analytic and satisfy the condition \2.1l . 
Now 

f^j |£«f + W)| 2 d*dy = fj |£f +1 «f +, (^,r ? )|d^drj. 

Here 

£f +i «f +i ^, I7 )^(£ M r i )(x, y ) V / }r T . 

Now 

+ Ej ?+l w 71 +F l p+l w 1 
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where the coefficients of the differential operator are analytic (smooth) functions of % and 



77 . Let Af be the unique polynomial which is the orthogonal projection of A p t into the 
space of polynomials of degree W in | and r\ with respect to the usual inner product in 



H\S). We define Bf +1 ,Cf +1 ,Df +1 ,E[ +1 and F ; p+1 in the same way. We then define 



p+i 



(£? +l y W =A 



p+i 



IB 



p+K 



-C 



"TJ7] 



D 



?p+i 



J?p+i 



Now let 7 be a side of the element £2, 
the mapping M„ + . Clearly 

j p+i 



P+1 



and let it be the image of the side % = under 



dx 



We now define 





i 


y dx 





P+U 



Here £, x (Q,rj) and rj x (0, 77) are the unique polynomials which are the orthogonal projec- 
tions of §c(0, 77) and 77^(0, 77) into the space of polynomials of degree W in % and 77 with 
respect to the usual inner product in H 2 (I). In the same way we can define (du{^ /dy) a 
on 7. Now let 7; be a side common to £2„ +1 and £2 P+1 and let it be the image of | = 
under the mapping Mj', + l and the image of | = 1 under the mapping M^ +1 



Let [w] denote the jump in w across 7/, where w is a smooth function on £1^ and 



ST 1 . We now define 



and 



/ du 
[dx 



( du 
\Jy 



1/2.M 



(5m 



p+i' 



1/2,7, 



<9.Y 



du?„ +1 
dy 



(0,T?) 



<9m„ 



p+i' 



(0,T?)- 



<5y 



(M) 



1/2,(0,1) 



(m: 



1/2,(0,1) 



Finally we consider a side of the polygonal domain £2 as shown in figure^ Let 7/ be a 
side of £2m + such that 7 C and such that 7 is the image of £ = under the mapping 
Mm +1 and which maps the master square 5 to £i£ + . Then we can define (du%f /dT) a 
and {du^ l /dN)" A in the same way. Finally we define 

2 

(0,17) 



p+i 



1/2,(0,1) 



and 



du \ ' 



1/2,7/ 



P+1' 



duf n 

dN 



(0,17; 



1/2,(0,1) 
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Now consider the sectoral domain £l k . Let us define the differential operator 
as in 1 3 1 . Then 

Shvfa, e k ) = a k w XkTk + 2p k w Tkdk + y k w BA + 8 k w Tk + e k w 6k + ^i k w, 

where the coefficients of £* are analytic functions of their arguments. Consider the ele- 
ment Clfj with 1 < j <M. Now the image of Cljj in (v kl fa) coordinates is the rectangle 

Q. k j. Clearly 



(£, k w(T k , e k )) 2 dx k de k = j & j (&ljw(v k , fa)) 2 dv k dfa 



Here 



£$jw(v k , fa) = £ k w(T k7 k )^Jj M k(v k ,fa), 



where J M k denotes the Jacobian of the transformation M k defined in J2.10I . Once more 
we can define a differential operator by replacing the coefficients of 2 k • by poly- 

nomials of degree W in v k and fa which are exponentially close approximation to them. 
Now the highest order terms of the differential operator £ k are given by SDt*, where 

Here y\ = T k and yi — 9 k . Let A k denote the 2x2 matrix such that A* > = 5^ •. Let // be a 
side of the element Q. k • such that // C T k , where F k is a side of the polygon fi. Let ji be 
the image of Yi in {yi ,y2j coordinates given by y\ = y\ (a), and y2 — yi(o)- Let f and n 
denote the unit tangent and normal vector at a point Pony/. We now define the conormal 
derivative 



dw 
dn 



J = n'A k V y w. 



Now the transformation M k defined in d2 . 1 Ot maps the rectangle Q. k j to -. Once more 
we can define (dw / dn)~ k \i^ by replacing the coefficients of the first order differential 
operator (dw/dn)^ k by polynomials of degree W in v k which are exponentially close 
approximations to them. We can now define || (dw/dn)^ ||^ 2 - as we have done before. 

The reader can now proceed directly to the stability theorem 13.21 stated in £|3 .31 and 
examine the proof later. 



3.2 Technical results 

Consider some Hf +1 € O p+l , as shown in fi gure|3] Then Glf + is a curvilinear quadrilat- 
eral whose sides are analytic arcs and the boundary <9£2j ,+1 is traversed in the clockwise 
direction. 
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Let 7 be a smooth curve and let N and T denote the unit outward normal and tangent 
vectors to 7 at a point P on 7. Let s be the arc length measured from a point on the curve 
in the clockwise direction. Then the second fundamental form is given by 

m^,n) = -^-T^ = ^-N^ = ^n, (3.8) 

OS OS 

where 

AT 

K = ± 

ds 

is the curvature of 7 at P. Clearly Trace(QS) = K. 

Now we need to use Theorem 3.1.1.2 of |4|. Let v be a smooth vector field defined 
on where v = (yi^vij ■ Consider the restriction of v to the boundary d£lf +l . Now 
d&f +1 = (U/=i 7')U(Uf=i Qi), where 7 are the sides of with end points deleted 

and Qi are the vertices of . We shall denote by vj the projection of v on the tangent 

vector T to <3iif + ' except at the vertices where this cannot be defined. Similarly by v^r 
we shall denote the component of v in the direction of N. Thus we have 

vn = v ■ N 

and 
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512/? 4 2 
^ f ( du 





du 


2 4 f 




dx r 


dx + 2£ / 




j=i J rj 



d / du 



ds 



<9m 



v'+i 



f du \ I du 
~ \dN~i) A \dfJ 



(Qj) 



(3.9) 



We shall say that a bounded open subset of R 2 with Lipschitz boundary F has a piecewise 
C 2 boundary if F = Fq \J Fi , where 

(a) To has zero measure (for the arc length measure ds) 

(b) Fi is open in F and each point x € F\ has a C 2 boundary as defined in 1.2.1.1 of |4|. 
Then Theorem 3.1.1.2 of |4| may be stated as follows: 

Let O be a bounded open subset of R 2 with Lipschitz boundary F. Assume in addition 
that r is piecewise C 2 . Then for all v S (H 2 (Q.)) 2 we have 



|div(v)| 2 dx- f V ^dx 
Jo r %dx s dx r 



r, \ds 



(v N v T )-2v T -v N \ds- l_ {(tr58)4 + 58(vr,v r )}ds. (3.10) 



To apply (13. lOi we define the vector field 



v = AV x m, 



where A is the matrix 

(A) r .s = a r 
We then observe that 



d ( d 
Mu= £ — [a r ^ ) Oliver). 



du\ du 

w A = L Nrar > s ^ =m ' N 



(3.11a) 
(3.11b) 



and 



du \ du 

— = £ T r a r , s — = (%v)-T. 



(3.11c) 
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Hence (13. lOi takes the form 



du \ 



if ki 



du 
~dN 



du \ 
df) 



dT J A ds \dNj A 

2\ 



ds 



ds. 



Now by Lemma 3.1.3.4 of |4| the following inequality holds for all u S H (£2) 



mo 2 L 



d 2 u 



dx r dx s 



2 

I 

r,s,k,l=l 



d u d u 



< Y. a r kd s i 3 — = — ^ — , 
■'-* dx s dXk ax r oxi 



a.e. in £2. Thus it follows that 



mo 2 L 

r,s=l 



d 2 U 



dx r dx s 



dx s dx r 



r,sM=\ 



a,-k 



d 2 u da s j du 



dxsdxk dx r dxi 



a.e. in £2. Integrating, we have 



mo 2 L 

r,s=l 



d 2 u 



dx r dx s 



dx< £ 

r..s=l 



5v r 5v s 
dx s dx r 



dx 



32R Z 



, 2 


(9m 


2 


5 2 M 




5x r 


I 

r..v= 1 


dx r dx s 



Ax 



where R is a common bound for all the C 1 norms of all the a,-. s . Hence 



(3.12) 



Next 



2 2 



r,5=l 



yi /• d 2 u z < ^ rdv r dv Sdx [ 512fl 4 2 
7 <5x r 5jc s ~ ,~-J dx s dx r " - i£ - 





(3m 


/ 


(9x r 



ds 



'v;yvr)ds = 



5m 



<9// 



5m 



/ 5m \ 



(2;)- 



dx. 



(3.13) 



(3.14) 



Then combining ( 13 . 1 2t — J3 . 1 4-1 we obtain the result. 

In a neighbourhood of the vertex we move to polar coordinates. We take a curvilinear 
rectangle £2* ■ which comprises part of the sectoral neighbourhood £2* of the vertex A^ 

and consider its image £2* • in (t*, 9k) variables as shown in figure^] 
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As in 1 3 ] we write the differential operator 9JI in modified polar coordinates, where 

r s= 1 r \ 5 / 



Now 



and 



X\ = Xj + e T * cos 0^ 



X2 = x\ +e T *' sin 9 k . 
Here A k = (jcj,.^). We would like to obtain an estimate for 



n 1 : 



rl\Wlu\ 2 dx: 



\Tl k u\ 2 dT k d9 k . 



Let us define the new differential operator 
2 d ( du 



J^, d ( du 

L ~ I a >- 



Here y\ — T k and yi = Q k - Let <9 A denote the matrix 



O k = 



cos k — sin 
sin k cos 0* 



(3.15) 



(3.16a) 



and A denote the matrix 

a i,l a i.2 

^2 1 ^2 2 

Then it can be easily shown that 



A k = 


(0*)'AO*. 






3 A 

n \ 






Q 2 r 5 ^ 






V 


A 
















Qi 


n 







(3.16b) 



Figure 4. Element £2*.. 
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Hence, since O k is an orthogonal matrix, we have that 

2 

£ ^ s T] r r] s >ii^l + nl). O.I?) 

r,s=l 

Moreover the following relations hold: 

(2t,iH=22i 2 + 0(e T *), (3.18a) 

(2t, 2 )%= 2*2.2 -5t,! + 0(e T *), (3.18b) 

(5* i2 ) % = -2a*,2 + 0(e^), (3.18c) 

(2t,i)r„(2t, 2 )% and (5| 2 K = 0(e T *), (3.18d) 
as Xk — > — 00 ■ Next let y be a curve given by 

x\ =xi(s), 
x 2 =x 2 {s), 

where s is the arc length along the curve y. Then the curvature K at a point P on the curve 
is given by 

dxi d 2 x 2 dx2 d 2 xi 



di di 2 di ds 2 
Let 7 be the image of the curve in (y i,y 2 ) coordinate given by 

yi =yi{a), 
yi=yi{(y), 

where a is the arc length along the curve y. Then it is easy to verify that 

i. (3.19) 



d^ 
da 



Now we can show that the curvature K of the curve y is given by 
da 

Hence 

\k\ < |K-|e T * + l <K, (3.20) 

where K is a uniform constant, for all the curves y s C £l k . 

We shall denote by t and n the unit tangent and outward normal vector at a point P on 
y, the boundary of Cl k j except at its vertices where these are not defined. 
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Lemma 3.2. Let u(y) G H 3 (£l k j). Then 



2 2 



2 



a 2 M 



dy r dy s 



dy 



< 



if / 3k \ f du 



<9m \ 



(6;) 



<9« 
d7 



<9m\ 

dn) 



da 



512 



* 4 E 



<9m 



5y r 



dy. 



(3.21) 



Now once more we use Theorem 3.1.1.2 of j4). Clearly Glf, for y > 2 is a bounded 

open subset of M. 2 with Lipschitz boundary T that is a piecewise C 2 . Thus r = ( IJ/Li 7i) 

U ( Uf=i Qi) where Ji are the sides of the open rectangle Q. 1 - • with the end points removed 

and 2i are its vertices. 
Now 

rl\Wlu\ 2 dx= [ e 4Tk \mu\ 2 dx k de k = [ \M k u\ 2 dy. 
in* ,. in* ; 



Here 



r f-j (9y r V r,S dy : 



du 



as defined in (I3.151 . Then for all w € (// (£2f •)) we have 



|div(vv)| 2 dy- V / ^l^idy 



4 

r 

>=l' / 55 



Here w„ and w, are the projections of w on the normal and tangent vectors n and f respec- 
tively. We define 



W = A VyU. 



Then 



du 



~k du 



[W 



(3.23a) 
(3.23b) 
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and 



(3.23c) 



So <I3 .221 takes the form 



l^yk l2j ^ <?W r <?W,s 



d f f du 



dn J i k 



da 



4 

L 



du\ 
dtJxk 



du 

dn j jt_ 



da 



4 

L 



du 



( du 



,dtJ+ l ) Ik \dnJ+ l 
Now using Lemma 3.1.3.4 of |4| we obtain 



du \ f du\ 
dil) Xk \d~nl) 



(Qj). (3.24) 



Mo 2 L 



d 2 u 



dy r dy s 



y dw r dw s 



2 

;„v.f./=l 



dy^dy, dy r dy/ 



and by J3. 1 8ai — J3. 1 8di there exists a constant R such that R is a common bound for the 

dw r dw s 



C l norms of all af ■. Hence 



u 2 2 
Mo y 

2 rh J % 



d 2 u 



dy r dy s 



Mo 2 ^ 



r=l" 





5m 


4, 





dr. 



(3.25) 



Thus combining ( 13.221 . ( 13 . 24i and ( 13.251 we get the result. 
We now need to write terms such as 



2p' 



Tfj 



du\ d / du \ 
dTj A d~ S \dN) A 



ds 



in (I3.2H where jj C = {(;ci,JC2) : = p} in terms of (yi,y2) coordinates. Let 7 be a 
smooth curve in £2jj = {(jci , JC2) : (x\,X2) 6 £2 and p<. < /x}, where p < /I, and let P be a 
point on 7 such that P in polar coordinates has the representation (p k , 9 k ) with p k = p. 
Now 



e?iV x u = 0*V y u, 



(3.26) 



where 0* is the matrix defined in (I3.16a> . and 



0% 



N = O k n. 



(3.27) 
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Hence 

e» ( j£) (P) = t'(O k )'AO%u(P) = t'A%u(P) = (P) (3.28a) 

using (13. 16al >. ( I3.26t and J3.27b . Here P is the image of the point P in {yi,yi) coordinates. 
Similarly, we have 



PROPOSITION 3.1. 

77i«s we can conclude that 



I ,du\ i,du\ AW da ( 3.29a) 

Jrj\dTj A d s \dNj A Jy } \dt) Ik da\dn) Ak 



and 



du\ fdu\ \ (p)= ((du\ (du 



In the same way we obtain the following results. 
PROPOSITION 3.2. 

Consider the boundary J common to Q k iM+l and £l k M . Then the following relations hold 
(figure^: 

4£M&)>>-{(£)*(&)>>- i3 « 



and 



Now let 7/ C 5 £2* , for some j <M and further suppose yi C where j £ Let ?z 
and f be the unit outward normal and tangent vectors, respectively, defined at every point 
of Yi ■ Then 
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X 



t 3 Q, 




n Q 2 



Figure 5. Elements Qr iM and £2* 



n 3 \tL 










N 2 


k 

T i,M+l J 4 




Y=y 2 [ 











- N 



M+r 



Here <T is the arc length measured from the point G (figure|6j where 

(t'A k n{a)) 2 



t(a)=t>A k t(a) 



n'A k n(a) 



and 



Hence 



h\a) 



t'A k n(a) 



n'A k n{a) 
du\ d f du 



dt ) da \dn 



da = If {a) 



Yl 



dud f du\ 
dt da \dn J 



da 



A* 



+ 



h k (a) d ( ( du 
dn 



iy t 2 da 
And so we can conclude that the following holds. 
PROPOSITION 3.3. 

d I du s 



da. 



4 do V dn) ^ k dCT 



^,.du d ( du\ ^ 
7l dt da \dn J 



1 f dh k (du* 2 



2 J % da \dn J 7k ' 2 



da 



h k (a) (du 



dnJ Xk 



Here g*(<7) and h k (a) are defined in J3.3 lbl and J3.3 lei . 



(3.31b) 



(3.31c) 



(3.32) 



Next let y m C dCl* , for some j > M such that y m C Tj where j G @. Let N and T be 
the unit normal and tangent vectors, respectively, defined at every point of y m . Then 



du 
df 



(ff) (s)+h(s) 



du \ 



(3.33a) 



Stability estimates for h-p spectral element methods 



415 




Figure 6. Arc length measured from the point G. 

where s is the arc length measured from the point G as shown in figure|6] Here 

(T'AN) 2 



g(s)=T t AT- 



N'AN 



and 



h(s) 



T'AN 



N'AN 

So we obtain the following result. 
PROPOSITION 3.4. 



P 

= P 



«9r; A d, V^A 

(9m d f du\ 
7 / {s) dTds{dN) A dS 

p 2 r dh (du\ 2 p 2 h (du\ 2 
Yj 7m ds~{dN) A S+ ~{dNj A 



(3.33b) 



(3.33c) 



(3.34) 



Now by J3.28bt we have that 



du \ " 

SvJ, 



And moreover by J3. 16al and 13.271 

*(G)=«*(S), 

and 

ft(G) =P(G). 
We can now prove the following estimate. 



(3.35a) 



(3.35b) 
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Lemma 3.3. Let u 1 ^ 1 G ff 3 (ilf + ). Then 

£ / /|Z)W M f +1 (^rj)| 2 d§dr, 

|ce|=2- ys - y 

■C^E IJ\DfD^uf +1 \ 2 ^dn\ 



<K /|£f +1 M f +1 | 2 d^dT ? +2p 2 £ 



'3 P+l 1 



d /<3«f +r 



/ ds\ dN 

/A \ 



Ep 2 



du, \ / oil. 



dN r+1 / I dT r+1 

/ A \ 



E / klP 2 



duf +l 
dN 



duf +1 
~df~ 



duf + 
dN r 

2^ 



duf + 
dT> 



ds 



(Qr) 



ds. 



(3.36) 



i/ere S is the unit square and £f + is the differential operator £ written in {£,,r\) coordi- 
nates. Here K and C are positive constants. 

Recall that 



2 2 

£m = - £ (flvW«nk + E b r (x)u Xr +c(x)u 

r,s= 1 / = 1 

= wiu+mu, 



(3.37) 



where 



Hence 



91m = E b r {x)u Xr +c(x)u. 

r=l 



p 2 f \Mu\ 2 dx<2p 2 f \£u\ 2 dx + 2p 2 f \mu\ 2 dx. 

Using Lemma l3~T1 we can conclude that there is a constant C such that the following 
estimate holds. 



9 9 9 
P Mo y 

2 /<->''+ 



«? 2 M f +1 



dx r dx s 



Ax 



cp 2 E 



r=l \ J "z 



dUf 



dx r 



dx + 



, |"f +1 | 2 dx 
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♦mm m-m mu 



J- 1 K\ / A \ / A \ I A 



dN \ dT 

/ A \ /A, 



As. (3.38) 



Writing the above in , r\ ) coordinates we obtain the result. 

In the same way we can prove the following estimate. 
Lemma 3.4. Let u k j <E H 3 (Q. k j). Then 

\a\=2- nl Ui J 

\\a\=l Jil iJ J Jil u J / 



< 



4 r 

r=VHr 



:*)>(*)> 



Heredtj = (\j/ k x (o!y,0!y +1 ) and fi ,C and K are positive constants. 

For 

£* M = e 2yi - £ (fl r ,v(x)M^) lr + Y, b r {x)u Xr +c(x)u 

\ r,s=l r=l J 

Here 

_ 2 _ 

m k u=Y,b k r (y)u yr +d k (y)u (3 .40) 

r=l 

and y = (yi ,^2) = (lib f° r some L Moreover the coefficients of 9^ satisfy 
^ = <9(e T *)forr= 1,2 
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and 

^ = 0(e 2T *) 

as — > — oo. 
Once more 




Using Lemma l3~2l we can conclude that there exists a constant C such that the following 
estimate holds. 




Rewriting (I3.4H in (v^, 0*) coordinates J3.39l > follows. 

We now need to obtain estimates for the spectral element functions in the H l norm 
which we do in the following theorem. 

Theorem 3.1. The following estimate holds: 

p h p m h l 

EEl«$il 2 +IIEll"^,fc)ll? ff +LH + (^)\\ls 

k=li=\ k=\j=2i=l ' '•' 1=1 

!p M h 
EEEll^/^Wfo 
fc=ly=2f=l ' W 

+ E E (IIMIIas + IIKlllo^ + IIKlfe) 

+ E E E (Ks+Kfe) 

ie9k=i-i ncda* f)r,, n(%)<oo 
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+ 1 L (ll[»]fe+IIK]||g, s +||[«fc]fe) 

*=ly s efl* 




(3.42) 



i/ere Cm = CM 4 f/;ere ex/sfs a vertex Aj such that Neumann boundary conditions are 
imposed on the adjoining sides T j and Tj+\ and Cm = C otherwise. C denotes a constant 
and n(Ys) the length of%. 

To prove the estimate 13.421 we shall use d3.5l >. To do so we have to define a corrected 
version of the spectral element functions so that it is conforming. 

Let {{u l l j (v k ,^ k )} i j<M,kA u Lj(^ r l)}iJ>M,k}k be a set of spectral element functions 
£ n M ' W . Here n M ' W is the set of spectral element functions such that u k ix — g k , a constant 
for all i, and u\ is a polynomial of degree W in each variable for j > 2. Then there is a 
set of spectral element functions 

{A*.(v*, <k)}ij<M, k , {A^(§ , i7)} iJ >M ) * e n M ' w 

such that the function (p(x\,X2) defined as 

(p(x u x 2 ) 




(rfj + tfj)(Vk(xuX2),<k(xi,X2)) if (xi,x 2 )€^j for j <M 
( uk ij + tfj)(Z(xi,x 2 ),T}(x 1 ,X2)) if (xi,x 2 )eSllj for j>M 



is a differentiable function of its arguments and <p G //q(£2). This can be shown as in 
Lemma 4.57 of Q. 
Moreover the estimate 

p h p m 4 

III^i| 2 +£ 1111^,^)11^, 

k=li=\ k=\j=2i=\ ' l 'j 

jfe=l ;=M+1/=1 
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< c I E E (ll«llo,s + ll M vJlo,r s ) 



K le®k=l-l y s CT,ndQ. k ,tl(%)<° 

E E (IIMfe + IIKlfe 

*=i»cn*/i($)<.= 



[«*Jfe) 



I E (IIMIIo.^ + IIKlll^ + IIWII^) 



E ll[«]ll§, yj - 



du 



dT 



-E E 

/e^%caop+inr; \ 



<9m 



«9r 



(3.43) 



0,% 



holds. 

We now explain the notation we have used in d3.43i . Let da denote an element of arc 
length in {Vk^k) coordinates. Then 



da. 



Moreover if y s is given by% = dQ£ +1 n^n +1 then 



du 
dT 



' 3 P+ 1 n P+ 1 ' 
a u m au„ 



dT dT 



ds. 



0,y s K 

Here d /dT denotes the tangential derivative in (jci ,X2) variables, i.e. 
5 m 



= T'V x u. 



The other terms in the right-hand side of (I3.43> are similarly defined. 
Now consider the bilinear form 



2 2 \ 

fl^WfcVj, + 52 M*)9'*r v + c( P v dx 

r,.v= 1 r= 1 / 



p IS. h. L 

jt=i j=i j=i ,J /=i 



Here 



2 2 \ 

52 flv(jc)<P% V ^ r + E b r(x)<Px r V + C(pv dx, 
r,s=l r=\ I 



B(<P,v)a = 

JA \r,.s = l 

where A is a domain contained in £2 and v 6 Hq (£2) 
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Now 



#0>v) n p+l = / „ +1 E a r,s( x )9x s V Xr + Y, b r{ x )9x r V + CfV dx 
1 Jn i \r,S=i r=\ J 

, , £(pvdx+ f , vds. 



Similarly if 1 <;'<Mwe have 



B((p,v) a * = A £ k (pvdT k d9 k + f ) vda. 



Moreover if j = 1 , 



B(9,v) rf = / c<pve 2T *d<r*d0* + / f ] vda 



since <p is a constant on j . 
Finally if j = M + 1 we obtain 



B(^,v) ot = / £<pvdx+ / I ^ ) vda 



Forby (E28B 



/ 



vdi. 



and ds = pda. Here P is any point on the circular arc B k p and P is its image in (t^, %) 
coordinates. Now 



P m 4 i 

Ar=l J=l 1=1 ' J /=! 



4-1 



p M h L 

fc=l./=li=l ' J Z=l ' 

( p M l k L \ 

\*=l/=li=l J /=1 ' / 

P M h i- _ Z. /■ 

jfc=ly-=l/=l-'l* f J /=!''"/ 



"I I /j^-^l 

*=ifccn*#(7»)<- * 



vda 
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+ E / 



du \ 



vd? +E E 



<9m\ 



vd(7 



EE E 



EE E 



(p 4 P M 4 

fc=l 1=1 ' Jfc=l ;=2i=l J 



+ LB(V fl ,v) or 
z=i ' 



(3.44) 



Now 



Here 



— ifr, or r, +1 cri°] 



O, otherwise j 
Now q = c(A k ), a constant, and c(xi ,x{) is an analytic function of x\ and X2. Hence 



£ k X^vdz k d9 k 



<2c k y lA^V^dTfcde, 

Jv 2 ^dX k dd k 



1/2 



for M large enough. And so we obtain 



L 



<e|A^|||v(xi,x 2 ) 



lo,n*,' 



where e is exponentially small in M. Now, let 2< j <M. Then 



£* M *; ; .vdT fe de* 



<||^«t-(T fc ,0 fc )|| o ^J|v(T,,0 fc )|| o ^ 



Finally 



Sc. 



nf +1 



(£wf +1 )vdx 



< ||£«f +1 (xi,X2)|| on p+i||v(xi,X2 



l 0,£l' 



P+l. 
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Now 

P m 4 

IIIllvfrbMljSa* <««||V(X1,X2)||U- 
ft=l ;=2i=l ' '•' 

Here ,Km = KM 2 if there is a vertex Aj such that Neumann boundary conditions are 
imposed on the adjoining sides Fj and and Km = A', otherwise. K denotes a constant. 
Hence 

P M h 

LLI>to.fc)ll?a* <^llv(^i^2)||?, Q . 0.45) 

fc=l ;'=2i=l ' ,J 

Now using the trace theorem for Sobolev spaces we obtain 
p m h 

<:= l./=2 /=! ' 'J 



And so we can conclude that 



p m 4 t 



(3.46) 



Using the Cauchy-Schwartz inequality in J3.44I and using ( I3.45> and ( I3.46> we can con- 
clude that 

{p m 4 „ p 4 

E E E ll£M;M)iioo* + E E e Kil 2 
fc=l ;=2(=1 ' 'J fc=li=l 



E E ... 



dcT+ £ 



* • / r» 



E E E 

/6^*=/-i r ,cr,ri5ii* 



du 

dnj~ Ak 



da 



£ / / |£ M f +1 (x 1 ,x 2 ) | 2 (kKk 2 
(9« \ 



y s cof 



dN 



ds 



E E 



i9m \ 2 



ds 



A*- 



dcr 



p 4 p m 4 

■eLLl\*il 2 +LLLll^AA)ii? ffi 



Jfc=l i=l 



k=l 7=2 1=1 
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P h 



P M 4 



p+ 1 



LLlN*i,*2)ii 



o,n* 



EEEiw^n 



l.nf 



Jfc=l;=2i=l 
L P if h f 

-Ellvr^ll^i + ELEL Vda 

/=! ' ' k=lj=2i=\ Jda! ij 



E y^+E E , 



v 2 ds 



Now v G //q (£2) and £ satisfies the inf-sup conditions (3.4). Hence using fOl i3~43l and 
(I3.46t we obtain 
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L 

E 



. , |£Mf +1 (x 1 ,x 2 )| 2 dxidx 2 



+ E (IIMIIo, ys +ll[^i]llo, rs +IIK2]llo, ys ) 



-E E ( ll«HS.»+ 



du 



dn 



dT 



0,7s, 
p 4 



E E y (§) d, +£ EE(i«?,i 2 +i^i 2 ) 



Here e is exponentially small in M. 

Using J3.43b and J3.45b once more we obtain the result. 
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We now define differential operators which are second order differential opera- 

tors with polynomial coefficients in and 0^ of degree W such that these coefficients are 
exponentially close approximation to the coefficients of (£*,•) as has been described in the 
beginning of this section. In the same way we define the differential operator [du/dn)~ k 
to be a first order differential operator with polynomial coefficients in and (j)^ such that 
these coefficients are exponentially close approximations to the coefficients of (du /dn)-r k . 
The other approximations are similarly defined. 

From the above, it is easy to conclude that 

Phi m \ 

EE l"ul 2 + EK;(v t ,0*)ll? fit 

k=li=\ \ j=2 ' V ) 

+ tK + \^il)\\l S <C M (y), (3.47) 

i=i 



where 
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+ E E E 
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1*110,54 + H H Xtllo,%. 



dn 
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Ell« 1 )X +1 (l,'7)ll^ 



+ E (l!Mllo, rt + l!Kri!^ + IIK 2 ril^) 



E E 



<9m 

df 



0,7s, 



E E 



du \ " 
dNj A 



0,rs 



Here Cm is as defined in Theorem l3.ll 
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3.3 The estimates 



We now define the quadratic form 

r M ' w (Hj(v k ,m,ikM +1 ^>n)}i) 
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(3.48) 



We can now state the main result of this section. 

Theorem 3.2. Let f M ' W ({ulj(v k ^ k )} iJ>k ,{uf +1 (^,T])}i) be as defined in dl48l . Then 
for M and W large enough the estimate 



p h 



M 



EE Kil 2 +Ell"U v ^* 



ZH +1 (^n)\\l, s 



k=li=\ 



j=2 



1=1 



C M ,wr MM (Rj(v fc) o A i}, a . K +1 (<^ , *7)},) 



(3.49) 



holds for all {{< / (v,,^)} lJ , it ,{«f +1 (^77)} / } S 7T M ' W 

//ere Cm,w = C maximum (M 4 ,(\nW) 2 ) if there is a vertex Aj such that Neu- 
mann boundary conditions are imposed on the adjoining sides Tj and Tj+i and 
Cm.w — C(lnW) 2 otherwise. C is a constant, independent of M andW. 

Adding a weighted combination of P.36> . ( I3.39> and P.47> and using the techniques 
and results of 1 3 1 the result follows. 
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Remark 2. The stability theorem 3.2 holds provided the coefficients of the differential 
operator 6 C 3 (£2) and the curves 0*, l and !//■*. f defined in (I2.6ab . (I2.6b> satisfy M.ll for 
' 1 3. 

For problems with mixed boundary conditions the factor multiplying the right-hand 
side of ( I3.49> grows rapidly with M. This creates difficulties in parallelizing the numeri- 
cal scheme. To overcome this we make the spectral element functions continuous at the 
vertices of the elements. Let 7ty' W denote the space of spectral element functions which 
are continuous at the vertices of their elements. We define 7Cr?' to be the space of spec- 
tral element functions which vanish at the vertices of their element. We now need to state 
a version of Theorem 13 .21 when the spectral element functions vanish at the vertices of 
their elements. 

To do so, we have to prove the following result. 

Lemma 3.5. Let ufj(^ , Tj) be a polynomial of degree W in B, and 7] separately, defined on 
the unit square S — (0, 1) x (0, 1), and which is zero at all the vertices of the square. Then 
there exists a positive constant C such that 

Mis < c(|«k(§ ,n)\ls + \4M Mh)- ( 3 - 50 ) 



Consider ,77) defined on (0, 1) x (0, 1). Now ^(0,0) = 0. Hence 
^°)=/ ^(^°)d§'. 



And so we can conclude that 



|«t-(S,o)| 2 <S 



Mi 

-af(«,o) 



Integrating the above with respect to £ we obtain 



\ul^M 2 ^<\ 



Mi 



<^(|«t;(^n)ll,5 + K/li,5) 
by the trace theorem for Sobolev spaces. Again 



(3.51) 
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Therefore 
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drj. 



Integrating the above with respect to | and tj we get 

JJ |^.(^r ? )| 2 d^d7 7 <2^ 1 |«u(^0)| 2 d^ 



Mi 
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Combining the above with I3.5H we obtain the required result. 

Clearly Lemma l3~5l appries equally well to any of the function elements iq ,-(Vfc, <j>k) for 
2 < j < M, 1 < i < Ik, 1 < k < p, although with a constant Q which depends on k. Taking 
the supremum over the constant Q (as given in ( I3.50t ) we conclude that 

< C(\4j(y*>ti)\ift + \ u U v k^k)\\a,\ (3.52) 

for all function elements with 1 < k < p, 1 < i < 4, 2 < j < M. Here C, of course, denotes 
a generic constant. We can now state the final result of this section. 

Theorem 3.3. Let {{w^(v^, , T7)}/} belong to the space of functions 

7Tq' W which are zero at the vertices of the elements on which they are defined. Then the 
following estimate holds: 

p M 4 

IID"y(vjt,fc)lE & +K (^)llv 

k=lj=2i=\ ' l <> 

<C(lnW) 2 r M - w ({«f, ; (v,,^)} I - M ,{ M f +1 (^T 7 )},) (3.53) 

for M and W large enough. 

In the above uf i(Vfc, fa) is taken to be identically zero for 1 < k < p and 1 < i < I k . 
Combining the estimates (13.501 and (13.521 with the earlier results (13.531 follows. 

4. Conclusion 

We can use the stability theorem l3~2l to formulate a numerical scheme to obtain an approx- 
imate solution to the elliptic boundary value problem i2.1\ as has been described in 18 19 1 . 
For problems with Dirichlet boundary conditions we choose our solution to be a non- 
conforming spectral element representation which minimizes a functional which is the 
sum of the squares of weighted squared norms of the residuals in the partial differential 
equation and fractional Sobolev norms of the residuals in the boundary conditions and a 
term which measures the sum of the jumps in the function and its derivatives in appropri- 
ate Sobolev norms at inter-element boundaries. In a sectoral neighbourhood of the corners 
these quantities are computed using modified polar coordinates and in the remaining part 
of the domain we use a global coordinate system. This method is faster than the h-p finite 
element method as there are no common boundary values to solve for 1 8 9 1 . 

For problems with mixed boundary conditions we have to make the spectral element 
functions continuous only at the vertices of the elements. As a result the Schur comple- 
ment matrix has a small dimension and an accurate inverse can be computed. Hence the 
numerical scheme has a computational complexity which is less for finite element meth- 
ods. 

Moreover, the construction of a pre-conditioner for the Schur complement matrix 
is very simple unlike the case for finite element methods. In fact, for problems in 
three dimensions the construction of pre-conditioners for the Schur complement matrix 
becomes quite complex for finite element methods |6 1. 

Though the ideas in these papers deal with problems in two dimensions, they generalize 
to three dimensions. We intend to study these problems both theoretically and computa- 
tionally in future work. 
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